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(54) Interleaver device and method for interleaving a data set 



(57) The interleaver device comprises a data proc- 
essor (16) for running an interleaver (I" 1 ) under the con- 
trol of driving means (18), input means (1 6) for inputting 
the data set to be interleaved and output means (18) for 
outputting the interleaved data set. 

The driving means (18) include mapping process- 
ing means (20) for performing a set of bijective elemen- 
tary functions (q> n ) and supplying a mapping of the 
interleaver to the data processor (16) for interleaving the 
data set according to this mapping and interleaver defi- 



nition means (24) for supplying said mapping process- 
ing means, with a definition of said interleaver (I 1 ), 
expressed as a compound function (<PkO...O(p-|) of ele- 
mentary functions (q) n ). each elementary functions com- 
ing from said set of bijective elementary functions fan), 
for said mapping processing means (20) to perform 
each of said functions, compounded according to the 
interleaver definition (I" 1 ), and so providing said map- 
ping to the data processor (16). 
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Description 

[0001 ] The invention concerns an interleaver device for interleaving a data set comprising a data processor for run- 
ning an interleaver under the control of driving means and an input means for inputting the data set to be interleaved 
5 and output means for outputting the interleaved data set. 

[0002] Such interleaver devices are specially useful in mobile telephones. 

[0003] An interleaver is usually applied on informations transmitted on a radio channel between two radio stations. 
[0004] In other technical fields, interleavers are also used for instance on data stored on a magnetic tape or a laser 
disk. 

10 [0005] Turbo decoders have proved to offer near Shannon limit error protection in the case of white additive gaus- 
sian noise. One key feature of the Turbo encoder is an interleaver. 
[0006] Interleaving data with a depth of N consists in: 

- writing N consecutive symbols of data into a buffer of size of N symbols, i th written symbol is written at position (i- 
15 1), and then 

- read them in another order given by some permutation I of {0, N-1}. i 1h read symbol is read from position l(i-1). 

[0007] In other words if i is the original position of a symbol in an input block, then l" 1 (i) is its position in the output 
block, where I" 1 denotes the inverse permutation of I. 
20 [0008] Deinterleaving data with a depth of N consists in: 

- writing N consecutive symbols of data into a buffer of size of N symbols, i th written symbol is written at position l(i- 
1). and 

- then and then read them in another order given by some permutation I on {0, .... N-1}. r read symbol is read from 
25 position i-1. 

[0009] To implement the interleaver and the deinterleaver, function is needed that provides the mapping I. With 
some simple permutation I, this mapping can be computed by a simple analytical formula that can be quickly evaluated 
on any known per se processing means. For instance if the interleaver is a rectangular interleaver with L lines and C 
30 columns, we have: 

• N = L • C 

• vi e {0,..., N-1} l" 1 (i) = (i div C) + (i mod C) • L (1) 

35 \ 

[001 0] In this formula (x div y) stands for the quotient of x divided by y in an Euclidian division and (x mod y) stands 
for the remainder in the same division. 

[0011] The way this formula is to be understood is quite simple: a rectangular interleaver consists of an array of L 
lines and C columns. The input data is written along the lines and read along the column. 
40 [001 2] If / is line number (from 0 to L-1 ) and c is a column number (from 0 to C-1 ) of some written symbol then writ- 
ing along the line comes to: 

c = i mod C ( 2 ) 
/ = i div C (3) 

that is to say: 



45 



50 



55 



i=/-C+C (4) 

/ and c are the co-ordinate of the (i+1) th written symbol. 
[001 3] Reading along the columns comes to: 

r 1 {i)-/ + c*L (5) 

[0014] Note that as a matter of fact, equations (4) and (5) are similar when column and line roles are inverted. By 
substituting equations (3) and (2) in equation (5), we get equation (1) that shows that the I mapping can be obtained by 
simple arithmetic computations in the case of a rectangular interleaver. 
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[001 5] However interleavers so simple as the rectangular interleaver do not always fit the needs of a forward error 
coding technique. In particular, in the case of the turbo encoding, the rectangular interleaver shows very poor perform- 
ance. The interleaver to be used in the turbo encoder must be more random than the regular rectangular interleaver, 
and however not completely random as it must still keep some good properties of spacing interleaved symbols. 
5 [001 6] The simplest way to implement that kind of interleaver is to use a table where the i th entry in the table is the 
value of l(i-1). The size of the table is directly derived from N, that is to say that at least s(N) bits are needed to imple- 
ment the table where: 

s(x) = xflog 2 (x)l (6) 



log 2 (x) denotes the logarithm in base two of x, and Txl denotes the ceiling of x, that is to say the least integer not 
15 less than x. When N grows bigger, then the amount of memory needed to store the table also grows bigger. 

[001 7] The aim of the invention is to provide an interleaver device which can run with very little memory to store the 
definition of the interleaver, even if the interleaver definition is very complex. 
[001 8] To this end, the invention relates to an interleaver device as defined in claim 1 . 
20 [001 9] According to particular embodiments, the interleaver device may include one or more of the features defined 
in subclaims. 

[0020] The invention proposes a method to build interleavers that are not so regular as the plain rectangular inter- 
leaver and nevertheless use far less memory than a completely tabulated interleaver. 

[0021] In the case of a mobile telephone, due to real time constraint, it might be necessary to tabulate the inter- 
25 leaver of the invention to be able to use it, since processing the values might take too much time. However the definition 
can be very concise, thus holding little space in ROM. The interleaver table can held in RAM and computed off-line, dur- 
ing connection establishment. After call completion, the RAM can be reused for another purpose. Thus, even if the inter- 
leaver needs to be tabulated into RAM for effective use, this is all the same favourable, because the RAM can be used 
for something else when the call type does not need the interleaver. 
30 [0022] Another advantage of having a concise definition from which a very large interleaver can be generated, is 
that a great variety of interleavers can be defined in the ROM memory, and only one of them selected at connection 
time. 

[0023] In a particular embodiment, the interleaver definition could be transmitted over the air, as a parameter of the 
connection. The interleaver definition would have to follow the mobile telephone if during a hand over of the site where 
35 the interleaving is performed in the network is changed. Therefore, the next site would need to compute or select the 
table fast enough during the hand over preparation time. 

[0024] The invention will be understood more clearly on reading the following description, given solely by way of 
example and made with reference to the drawings, in which : 

40 - Fig. 1 is a block diagram of an interleaver device according to the invention; 

Fig. 2 is a flow-chart showing a Mod-threshold term splitting algorithm; 

Fig. 3 is an illustration of a term splitting function applied on a data set; 

Fig. 4 is an illustration of a Div-mod factor split function applied on a data set; 

Fig. 5 is an illustration of a mod-mod factor split function applied on a data set; 
45 - Fig 6 is an illustration of a rectangular interleaver applied to a data set; 

Fig. 7 is an illustration of a truncated rectangular interleaver applied to a data set ; and 

Fig. 8 is an illustration of a puncturing of an interleaver. 

[0025] The interleaver according to the invention can be implemented on any data processing means adapted to 
so run under the driving of a software using the splitting of the interleaver as defined hereafter. 
[0026] For example, the interleaver can be implemented in a mobile telephone. 

[0027] Generally, and as shown on figure 1 , the interleaver device includes a processing unit 10. It also includes 
means 12 for inputting a data set to be interleaved and means 14 for outputting the interleaved data. 
[0028] The processing unit 10 includes a data processor 16 connected to the inputting means 12 for receiving a 
55 data set to be interleaved. It is also connected to outputting means 1 4 for supplying the interleaved data. 

[0029] The data processor 16 is adapted for running an interleaved I" 1 , the mapping of which is supplied by driving 
means 18. The driving means 18 comprise mapping processing means 20 which are adapted to compute a mapping 
of the interleaver from a set of elementary bijective functions <p n which can be derived from a Read only memory 22. 
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The computation is performed, according to an interleaver definition I" 1 received from interleaver definition means 24. 
[0030] In the present embodiment, the elementary functions are stored in the Read only memory 22. Typically 
these elementary functions may be not completely defined in ROM 22, and zero one or more additional parameters may 
need to be supplied for them to be operated on their argument. 

5 [0031] Interleaver definition means 24, typically consisting of a memory, are provided to supply the definition of the 
interleaver I" 1 to mapping processing means 20. This definition is based on the elementary functions stored in memory 
22. In particular, the interleaver I* 1 is defined as a compound of the elementary functions <p n stored in ROM 22, along 
with the suitable parameters if any. The mapping processing means 20 drive the data processor 10 for it to interleave 
data according to the mapping which is computed by the mapping processing means 20. This mapping is determined 

10 by performing each of the functions <p n used to define the interleaver compounded according to the interleaver definition 

r 1 . 

[0032] Typically the mapping processing means 20 compute the mapping from the definitions from ROM 22 and 
interleaver definition means 24 only once during connection establishment, and then keep it in a look-up table. In 
another embodiment the mapping processing means 24 compute from the definition of I" 1 the image l" 1 (i) of the posi- 
15 tion i of a data symbol to be interleaved on the fly at every symbol. 

[0033] The implementation of the interleaver can be easily carried out by a man skilled in the art in view of the fol- 
lowing explanations. 

[0034] The definition of interieavers by the inventive method is recursive. The interleaver is defined by a permuta- 
tion r 1 on the set {0, 1 , .... N-1) of the initial position indexes of the data to be interleaved. This set of position indexes, 
20 on which operates the interleaver I is split into a sum or a product of two smaller sets. Each of these two smaller sets 
can be then : 

either permuted by the underlying index permutation of a smaller interleaver, 
either split itself again in a sum or a product of two other smaller sets, or 
25 - merged with a smaller set that was produced during a previous split of set. 

[0035] The invention allows the use of interieavers external to the invention in place of the smaller interieavers. 
These external interieavers can be tabulated, or computed by a processing method different from that of the invention. 
[0036] In the description, the following notations are used. 

30 

FINITE ELEMENTARY SET OF INTEGER 

[0037] For any non null positive integer x, H denotes the set: 

35 — — 

[x] = {i/ i integer and 0 £ i £ x-1} 



40 

ELLIPSIS NOTATION 

[0038] The ellipsis (...) notation is used in the middle of a sequence instead of the full sequence. This notation does 
not make any assumption on the number of literals in the formal sequence using it. These convention are quite obvious 
45 and usual. Below, some instances are given: 



so 
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10 



15 



• Aix... xA p Shall be interpreted just as Ai if p = 1 

• Ai#. . . 4*A P shall be interpreted just as Ai if p = 1 

• A^. . . i^Ap shall be interpreted just as the empty set 0 if p = 0 

• {x f x+1 y-1,y} where x and y are integers, shall be interpreted just as {x} if 

x=y 

• {x,x+1, .... y-1,y} where x and y are integers, shall be interpreted the empty set 

0 if x>y 



PRODUCT OF SETS 

20 [0039] If A and B are two sets, then their product AxB is the set of couples (a,b) such that a is in A and b in B. In 
other words: 

AxB = {(a,b)/a e Aandb e B} 

25 [0040] The concept can be generalised to any finite number of sets: 

A 1 x...xA p -{(a 1 a p )/via, e Aj) 

[0041] Furthermore, for the sake of simplicity, Ax(BxC), or (AxB)xC will just be denoted AxBxC, since there are 
30 obvious bijections between these three sets, that is to say there is little difference between ((a,b),c), (a,(b,c)) and (a.b.c). 
[0042] One property of the product of sets is that if Ai , A 2 , ... Ap are finite sets and their respective numbers of ele- 
ments are denoted l^l, |A 2 |, ... |Ap|, then A^... xAp is also finite, and its number of elements lA^... xAp| is such that: 

i= P 

35 lA^xApl-JlN 

1*1 



UNION OF SETS 

40 

[0043] The union (noted of two sets is the set of elements that are in at least one of these two sets. That is to 
say AuB = {x/xeAorXeB}. 

DISJOINT UNION OF SETS (ALSO CALLED HERE SUM OF SETS) 

45 

[0044] If A and B are two sets, then their sum A * B is the set of couples (t,x) such that t is equal to 1 or 2 respec- 
tively when x is A or in B. Consequently, 



50 



A*B = ({1}xA)u({2}xB). 



[0045] 1 and 2 are the usual integers and are used to distinguish the elements of A and of B in the union, so that 
the union is "disjoint". 
55 [0046] The concept can be generalised to any finite number of sets: 
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A^... tt>Ap = ({1}xAi)vj ... vj ({p}xAp) 



[0047] Where 1, 2, ...p are the usual whole numbers. 

[0048] In the following for (t,x) in A 1 * ... & Ap, t is called the tag and x is called the value. 

[0049] Furthermore, for the sake of simplicity, Av* (B^C), and (A^B) * C will just be denoted AtbBtbC, since 
there are obvious bijections between these three sets, that is to say there is little difference between the tag sets {1, 
(2,1), (2,2)} , {1, (1,1), (1.2)},and {1.2,3}. The main use of the tag is its ability to distinguish the elements according to 
their set of origin. 

[0050] Note that the disjoint union is generally defined in literature such that 

A^B = B^A p 

that is to say there is no underlying order on the set of tags, and therefore neither are the terms ordered. However the 
definition is slightly modified, because this property is not desirable for this application. In the following, the expression 
"sum of sets" instead of "disjoint union of sets" will be used. 

[0051] Finally, one property of the disjoint union of sets is that if A 1p Ag, ... Ap are finite sets and their respective 
number of elements are denoted \^\, |A 2 |, ...|Ap|, then A^.-.^Ap is also finite, and its number of elements 
lA^.^^Apl is such that: 

i=P 

lA^.^Ap^S^I 
1=1 



DISTRIBUTION OF PRODUCT OF SETS OVER SUM OF SETS. 
[0052] In the following it is considered that: 

Ax(B^bC) is the same as (AxB)^(AxC) 



[0053] As a matter of fact, there is an obvious bijection mapping the element (a,(t,v)) of Ax(B * C) on the element 
(t,(a,v))of 

(AxB)^(AxC). 

This bijection just consists in placing the tag in first position. 
[0054] Similarly it is considered that (A * B) xC is the same as 



[0055] The interleaver I is built as follows. It enables a concise definition of the interleaver using several bijective 
elementary functions. 

[0056] The general idea to is to build the bijection I" 1 as the compound of several bijections q> 1p 92 <pr- 

Thatistosay I" =<p k o....9 2 -9i 
[0057] This way we have: 
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vn G {1....k}<p n S (n - 1) ->S (n) 



s (0) = [l 



S (k) = [N]. 



15 [0058] This way we have: 



IS 



20 



S (0) cp^ S (D <p| 
x i— * cpi (x) 



— > ® 

<P£ S (k) 

T 1 (X)= q>k(q)k-l(...92(9i(x))...)) 



25 and I" 1 is so built as a compound or sequence of k simpler steps, each step consisting of a bijective function <p n . 

One more aspect is that each intermediate set can be defined as a sum of simpler sets Tf n ) that are products of 
sets of the H form: 



30 



35 



40 



45 



Vn €{1,...,k}S< n) =T < r ) *...*Tq n 
where : 



(n) 

In 



Vn e{1 k}Vie{1, q„} T, (n) = [n£? 



N|,2 



x...x 



N 



M 

i.Pl.n 



and with (for all the q> n are bijective), all the intermediate sets S (n) having N elements : 

z n 

i=1 V j=1 



vne{0 k},|S (n) |=X(n NlJJ-N 



and also q 0 = 1 Pi,o = 1. N 1 "', =N ,q k = 1, p 1<k = 1 and =N, because 

s (0) = s W = [n] 



so [0059] Note that for the present application, it is needed that each <p n be easily implementable by an algorithm that 
can be run on existing processing means. 

[0060] Several types of <p n elementary functions are defined. The different tp n elementary functions form a set of ele- 
mentary functions which are stored in memory 22, so that their algorithm can be run by the mapping processing means 
20 according to the interleaver definition I" 1 . 
55 [0061] Here is an introductory list of the functions that will be presented in more details later: 

A. Functions to simplify the definitions of the other functions by putting in first position the terms or factors on which 
the next function is to operate: 
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■ Term permutation 

■ Factor permutation 

B. Functions to decompose S*"* 1 ) into an with more terms or factors: 

■ Elementary term split 

■ Elementary factor split 

C. Functions ensuring the equivalence between Ax(B^G) and 

(AxB)^(AxC) 



■ Factorisation 

■ Distribution 

20 D. Function to interleave a part of S (n " 1) 

■ Elementary interleaving: embedding of an interleaver external to the invention 

■ Affine-mod: a linear relation followed by a modulus, and operating on a term of S< n " 1 ) as if it was a vector space 
whose dimension is the number of factors in this term, 



E. Functions to simplify S< n ~ 1) into an S (n > with fewer terms or factors: 

■ Elementary factor merge 

■ Elementary term merge 



[0062] The different functions are as follows. 
TERM PERMUTATION % 
35 [0063] In that case: 

" Q n = n-1 

a a permutation of {1,..., q n }, and 
- vie{1 q n }T^=Tr 

40 



9 n is defined as: 



45 



cp n : S (n - 1) S (n) 

(t,v) h* (o(t),v) 



so [0064] In the definition above t is the tag of (t,v). and v the value. 
Example: 



55 
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S (n> = (U * Bxgl * [2]x[3] 



<* = (3,2,1) 

ro And for instance we have the following mappings: 



15 



<p n : 



(1.(1,0) 
(3.(2)) 
(2.(3,0)) 



(3.(1,1)) 
(1.(2)) 
(2.(3,0)) 



20 



25 



30 



35 



40 



FACTOR PERMUTATION q> a The factor permutation <p n function permutes the factors in the product constituting the first 
term of S (n " 1) . In fact the same principle could be used to permutes the factors of any term, but this is not needed for 
the definition thanks to the term permutation function defined above. In that case we have: 



[0065] 



qn = ( Vi 

Pl,n "P Ln-1 

ct a permutation of {1 pj n .^} 



vi*{2 q n }Tf =T[ 



*M1 P w }Nf3.Nf;i 



9 n is defined as: 



x = ( t, ( x 1l x 2 ...x pM)) 



S (n) 

f( 1 '( X o(1)' X a(2)- X a(p 1 , n . 1 ))) ift=1 

Ix otherwise 



45 



50 



Example: 



S (n1) = H]x[3]x[4] * g) 
S (M, = [3]x[2]x(4] * g] 



55 



a =(2,1 ,3) 

And for instance we have the following mappings: 
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cp„: (1.(1.2,3)) w (1.(2,1.3)) 

(2,(2)) h> (2.(2)) 

(1.(0,1,2)) h> (1.(1.0,2)) 



10 

[0066] For the term or factor permutation, <p n is made at null processing cost. The only added value of <Pn 's that it 
makes simpler the definitions that will follow by alleviating the notation. The following definitions are applicable to any 
terms or factors in a sum or product of sets, respectively, instead, thanks to the factor permutation (p n we can without 
any loss of generality make them only for the first or the first and some subsequent terms or factors according to the 

15 case. 

[0067] In consequence the permutation cr is typically the compound of zero, one or two transpositions. A transposi- 
tion is a permutation that swaps two elements of a set, and let the other unchanged. Thus such a permutation cr puts in 
the first position the terms or factors on which the next function <p n is to operate. In other words the term or factor per- 
mutation function specifies on which term(s) or factor(s) is operating the next function in the <Pk°...°<j>2°9i compound. 

20 

ELEMENTARY TERM SPLIT <D n 

[0068] The elementary term split function q> n is splitting the first factor in the first term of S (n ' 1) into a & sum of two 
sets. In fact, the same principle could be used to split any factor of any term, but this is not needed for the definition; 
25 thanks to the term and factor permutation functions defined above. An example is give n on Figure 3. 

[0069] S<") is such that from S< n " 1) to S< n ), nothing is changed except the first factor IA+BI of the first term T A < n ' 1) that 
is split by a bijection f into the * sum of 0 and [E5| on which are then distributed the subsequent factors, if any, of the 
first term of 

[0070] In that case : 

30 

q n = q n _.,+1 , i.e. one more term 

• vi g {2,...,q n _i } T["{ = Tf n l) , i.e. the subsequent terms unchanged 

• A and B are two positive and non null integers, and ] = A+B , = A and = B 

• P 1in =P 2 ,n=Pi t n-i andvj G {2 P^IN^-nSJ-N^j 11 

35 • a bijection 

* |A+a -> 0*11 

40 x h-> (tKx).Vf(x)) 



45 where t^x) and v,(x) respectively denote the tag and the value of f(x) in 0 * H . 

[0071] Note that the conditions of 3 rd and 4 th bullets are such that : 

T ( r ,) = IaUx r 



55 
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w 



15 



20 



25 



<p n is built as: 

9n: S<"- 1 > 



with 
x = (t,v) 
and 



T ( i n> = 0 
And T = 



T 
T 



V= jX^Xj-.X^, 



N (n-D 



X...X 



N 1.Pl.n-1 



s (n) 

LtKxi^VfCXf). x 2 ...x piM)) ift=l 
l(t+1,v) otherwise 



[0072] The f bijections used by the elementary term split function can be defined as follows. 



30 



f : | A+B | -> [ a]^[b] main functions 



35 



40 



[0073] There are many possible functions that all are computable on any current processor. 

[0074] As an example of f, the mod-threshold term spli t func tion will be described hereafter: 

[0075] In this function, the modulus of an element x of |A+B[ by some constant C is computed. Then it is compared 

to a threshold T According to the result of this comparison, it is decided to map x either to E) or to [@. Constants C and 

T are such that: 

0<T<C 



[0076] In order to understand more easily how the mod-threshold function works, a simplified algorithm will be first 
presented. The simplified algorithm cannot make the mapping in a random way, but only in a sequential way, that is to 
45 say it can map element x, only after having mapped elements 0, 1, ...,x-1 before. 

[0077] In the simplified algorithm there are two counters n A and n B that hold the number of elements that have 
already been respectively mapped to [A] and H. 



50 
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n A := 0; 
n B 0; 

for x:= 0 to A+B-1 do 
begin 

if (n A = A) then 
begin 

remark [aJ is full, then map x to [b] 

'5 map x to (2,n B ); 

n B :- ne + 1 ; 
end 

20 else if (ne - B) then 

begin 

remark \b\\s full, then map x to [a] 

mapxto(1,nA); 
n A :=n A + 1; 
end 
else 
begin 

35 remark neither [a) or \b\ are yet full, then map x accord- 

ing to threshold of modulus 
c := x mod C; 
if c<T then 
begin 
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mapxto(1,n A ); 
n A := n A + 1; 



end 



else 



begin 

map x to (2,n B ); 
riB := Ob + 1 ; 

end 



end 



end 



[0078] Now, an algorithm that can work with random values of x in input (x not given in a sequential way) is pre- 
sented. 

[0079] Three constants x^, and tM are derived from A, B, C and T according to the following definitions and for- 
mulas: x M is the greatest x while neither 0 nor |@ are yet full in the simplified algorithm above. 

t M is the tag of the set (0 or [B|) not yet full after x M has been reached by x in the simplified algorithm 
v M is the greatest value of the set not yet full while both sets are not yet full in the simplified algorithm. 
x Mp tM and v M can be alternatively defined by the following formulas: 



A 

x M = max 



jx € [A+Bl / 



f/-T + c<A ifc<T 
and 

V /-(C-T) + c-T < B otherwise 



, where / = x div C, and c = x mod C 



t M -Q 



if ((x M +1)modC)>T 
otherwise 



ift M = 1 thenv M ^ max S v /3x <; x Ml 



^xdivC 
and 

Vc -x mod C 



'c<T 
and 
ly = IT + c 



ift M = 2thenv M ^ max"* 



/3x 



^ Xm, 



^xdivC 
and 

\c -x mod C 



and 

[y = / (C-T) + c-T 
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[0080] Then the algorithm that maps x onto (t.v) is the following: 



(t.v) 

defined by the algorithm below 

if (x ^ x M ) then 
begin 

remark Neither [A] nor [b] are yet full, then discrimination is based 

15 on threshold 

if (x mod C) < T then 
begin 

t:=1; 

v := (x div C) T + (x mod C); 

end 
else 
begin 

t:=2; 

v := (x dlv C) (C-T) + ((x mod C)-T); 

end 

end 
else 
begin 

remark [A] or [|] is already full, then no more discrimination, the 



remainder of |A+B[ is mapped to the set not yet full, that is 
45 to say that with tag t M 

t :=t M ; 

v := x - x M + v M ; 

end 



[0081] The Mod-threshold term splitting algorithm is also shown on Figure 2. The result of the algorithm is illus- 
55 trated on figure 3 in the case of A=5, B=10, C=4 and T«2. 

[0082] For some particular value to T and C the definition can be simplified. For instance when C = A+B and 
T = A , f is defined as: 
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f: 



A+B 



10 



[A]«© 
(t.v) 

with if x < A 

then t = 1 and v = x 
else t = 2 and v = x-A 



15 [0083] Another simplification occurs when given the values of D=A+B, C and T , A and B are chosen such that the 
comparison of x to x M is not necessary, and only the comparison of the modulus to the threshold needs to be done. That 
is to say the mod-threshold term split function can be more simply written as: 



20 



25 



30 



f: 



[MBJ 
x 



0*H 

(t,v) 

with if c < T 

then t = 1 and v = IT+c 
else t = 2 and v = / (C-T)+c-T 



where" 



= xmod C 
x div C 



[0084] Given D = A+B, C and T , the value of A and B for which this simplification occurs can be computed as: 



35 



40 



A=|{xe[Dj/ (xmodC)<T}| 
B= |{xe[Dj/ (xmodC£T}| 



(7) 



where |.j denotes the number of elements in a set. 
45 [0085] This can also be more simply computed as: 



50 



55 



A = 1 + / A T + c A 



where 4 



f I A = x A div C, 
c A = x A mod C, and 

x A = max{x e \5\ I (x mcxl C)<T } 

{/b = x B div C, 
c B = x B mod C, and 
x B = maxjx e {5\ I (x mod C)zT ) 
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[0086] Finally x A and Xb can be more simply computed by the following algorithm: 

c := (D-1) mod C; 
if c < T then 
begin 

x A :=(D-1); 

x B := (D-1)-(c+1); 

end 
else 
begin 

x B := (D-1); 

x A := (D-1)^(c+1)+T; 

end 



25 [0087] For instance D=1 0 C=5 T=2 will yield x B = 9, x A = 6, / A = 1 , c A = 1 , / B = 1 , c B - 4, A = 4 and B a 6. 

[0088] Selecting A and B such that the algorithm is simplified can result in less processing power requirement to 
prepare to interleaver table, or to compute the interleaver real time, or less Integrated Circuit surface, if this real time 
term splitting procedure is performed by an ASIC. 

30 ELEMENTARY FACTOR SPLIT g> n 

[0089] The elementary factor split (p n function is splitting the first factor in the first term of S< n " 1 > into a x product of 
two sets. In fact the same principle could be used to split any factor of any term, but this is not needed for the definition 
thanks to the term and factor permutation functions defined above. 
35 [0090] S< n > is such that from S< n_1 > to nothing is changed except the first factor Ia^bI of the first term that is 
split by a bijection g into the x product of 0 and [|. 
[0091] In that case: 

q n = q |v1 . i.e. same number of terms 
40 * vi e{2,...,q n . 1 } T| n) = T[ n1) , i.e. the subsequent terms unchanged 

P i n - P i n-i +1 « i e - one more * actor in f irst t er m 

• Pn-l}N^i=N^ 

• A and B defined by A = Nip) , B=N£0 , and thus A • B = 
a bijection 

45 

W (gi(x). g 2 (x)) 

50 

cp n is built as: 

55 



AB 
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; S (M) -> S (n) 



10 



15 



X 

with 
x = (t,v) 
and 

v= ( x vVV) 



l( 1.(9^), ga^), x 2 ... , x^^ if t=1 
lx otherwise 



[0092] The g bijections used by the elementary factor split function can be defined as follows: 

20 



fll |Aj3|-> [ A]x[b] main factor split functions 



25 

[0093] There are many possible functions that all are computable on any current processor. 

[0094] Here are some examples: 

[0095] In the following for all x and y integer such that y>0: 

30 - xdivy stands for the quotient of Euclidian division of x by y, 

x mod y stands for the remainder of Euclidian division of x by y. 

[0096] Note that vx vy y>0, (x mod y) e [y], even for x ^ 0 
[0097] The DIV-MOD factor split function: 



35 



40 



x j-» (x div B, x mod B) 



[0098] This function is clearly bijective and: 

45 g" 1 ((a,b))=a«B+b 

[0099] The div-mod factor split function is illustrated on Figure 4. 
[0100] The MOD-MOD factor split function: 

This function can be used only when the greatest common divider of A and B is 1 (that is to say A a B = 1). 

50 

g: [a! -> ExGH 

x k> (x mod A, x mod B) 

[0101] The fact that the greatest common divider of A and B is 1 ensures the existence of two integers A' and B' 



55 
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such that B • B'+A • A' = 1 . 

[01 02] This allows to build easily the inverse of g as: g- 1 ((a,b) ) = b + { ((a-b) • B') mod A) • B = (a • B • B' + b • A • A') 
mod (A • B) 

[0103] The mod-mod factor split function is illustrated on Figure 5. 

5 

FACTORISATION Vq 

[0104] The factorisation <p n function groups several equal terms of a sum of sets, into a product of sets. This uses 
the obvious bijection between: 

10 

B*_*B And \K\xB 



15 

A times 



20 [0105] This bijection simply maps the tag-value couple (t,v) on the value only list (t-1 ,v), making of the tag t an ele- 
ment M of 0. 

[0106] The factorisation q> n thus factorises the A first terms for S^" 1 ). In fact the same principle could be used to 
factorise any terms in S (n " 1) provided that they are equal, but this is not needed for the definition thanks to the term per- 
mutation functions defined above. 
25 [01 07] S< n ) and S< M ) are such that: 

• A is positive non null integer such that q n » q fr1 -A+1 , because A terms of S (n " 1) are merged in one, 

• vie{1 A} m T^ , that is to say the A first terms of S^' 1 ) are identical, 



30 



35 



40 



45 



vie {2,..., q n } T[ n * = .that is to say the terms following the A first terms of S*"" 1 ) are unaffected 
<p n is built as: 



S (rv1) ^ g<n) 



X 

with 
x = (t,v) 
and 

V = ( X 1' X 2- X J 



h> ( y1 l( Xj V t x pini)) ift=l 

l(t+A-1,v) otherwise 



DISTRIBUTION q> n 

[01 08] Trie distribution <p n function does the inverse operation of the factorisation <p*n that was described in the pre- 
vious section. That is to say the first factor of the first term in S (rv1} becomes a tag for a sum. 
55 [0109] In fact the same principle could be used on any factor of any term of S (n1 ) provided that the term containing 
that factor has at least two factors, but this is not needed for the definition thanks to the term permutation functions 
defined above. 

[01 1 0] S<") and S^ 1 > are such that: 
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10 



15 



20 



A is a positive non null integer such that N^ 1 * = A 

q n = q + A-1 , that is to say 1 term of S (n " 1) is distributed over A terms of S {n) 
Pi ( n-i > 1 ( f irst term of s(n ' 1) nas at ,east two factors) 

v i g {1 A} p | = p 1 n _., -1 , the A first terms of $( n ) have one factor fewer - the one that was distributed - than the 

first term of 

vie{1,...,A}vje{1 p 1 n -1 J N[J } = N^j , the Afirst terms of S< n > are all identical and equal to thefirst term of S< n ' 

1 ) without its first factor' 

vie {2 q n .^} = t|"' 1) p the q^-l subsequent terms of S (n * 1) are unaffected 

<p n is built as: 



SOi-D 
x 

with 
x = (t,v) 
and 



S (n) 



W J(V 1 '( X 2' X 3' 
l(t+A-1,v) 



Pl,n.l) ) 



if t=1 

otherwise 



25 



30 



35 



ELEMENTARY I NTERLEAVE ¥n 

[01 1 1 ] The elementary interleaving (p n function permutes the elements of the first factor in the first term of S< n " 1) . In 
fact the same principle could be used to permute the elements of interleave any factor of any term, but this is not 
needed for the definition thanks to the term and factor permutation functions defined above. 
[0112] S (n) and S (n1) are of course equal. 
[0113] In that case: 



s (n) = s v 

A an integer such that A = 
a bijection i : El -> G9 <P n ' s built as: 



5 <n-1) 



40 



45 



50 



55 



cpn: 

X 

with 
x = (t,v) 
and 



$ <n) 

.'(LtUX,). V ■Vn-l)) ift=1 

x otherwise 



[01 1 4] Here are some examples of i bijection used by the elementary permutation function. 
[0115] The tabulated i function: 
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[0116] In that case the i is implemented by a table in the memory of the processing unit. There are A entries in 
that table and » (x) is written in the (x+1) th entry. 

[01 1 7] The table therefore requires s(A) bits (where s is defined in equation (6)) 
5 AFFINE-MOD <o n 

[01 18] The affine-mod <p n function acts in the first terms of S fn ' 1) as if it was a vector space of dimension p 1 n _i . In 
fact the same principle could be used on any term, but this is not needed for the definition thanks to the term permuta- 
tion functions defined above. 
10 [01 19] S (n) and S (n " 1 > are of course equal. 
[0120] In that case: 



15 



20 



25 



s (n) = s (n " 1) 

m = Pl,n1 

an m line m column matrix U=[u , j] such that 

vie{1 m}vje{1 m} u g is integer (possibly negative) 

a m line vector V=[v ,] such that: 

vie{1,... ,m}Vj is integer (possibly negative) 

<P n is buiit as: 

cp n : S (n - 1 > -> S (n) 



30 



35 



40 



(t,x) 



((1,y)ift = 1 
l(t,x) otherwise 

with for t=1 

X — ^ X|, T X2- . 

y = ( y 1 .y 2 - 

such that: 









"xr 






=U- 




45 


-Z m _ 




_Xm- 



» x m) 
■Vm) 



and 



Vie{1 f .. M m}yi =z, mod Nf 



(n) 



50 



[01 21 ] Furthermore U and V are such that the tp n obtained be bijective. 
55 [01 22] One sufficient condition to be fulfilled by U and V for the affine-mod <p n to be bijective is the following one: 
U is the product of two matrices and U^. 

- The elements of these matrices are denoted this way: vk e {1 ,2} U (k> = [uff] 
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U (1) is a sub-diagonal matrix with a diagonal of 1 or -1 , that is to say: 



V (ij)e {1 m} 2 j>i=>u ( J=0 

and 

Vi€{1 m} ( u?/=1 oru ( J ) = *1 



- U< 2 > is a diagonal matrix with for all i , i* diagonal element primary with N M , 
[0123] That is to say: 



15 



V(i.j)e {1 m} 2 H=>u ( J=0 

and 



V ie {l.-.-.m}, and NVj greatest common divider is 1 



,{n) 



20 

- Product is done in this order: U = U {2) • U (1) 

Note: U (1) can be also super-diagonal (elements under diagonal are null) instead of sub-diagonal, this comes up to the 
25, same thanks to the factor permutation <p n function. 

ANNEXE ON AFFINE-MQD FUNCTIONS 

[0124] Here is given a proof of the sufficient condition that was given to build the affine-mod q> n functions. Let us 
30 recall the definition of the affine-mod function where we have omitted the n index that brings no information here, and 
we have replaced it by a U index, meaning that this is the q> that is generated by some matrix U. For alleviating the nota- 
tion we have also omitted the term index (specifying that we are working on the first term of S (n ' 1) ). 



35 



40 



45 
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( X V X 2" X p) 



[Ni|x|N 2 jx...x|N p 

( y v y 2 ...,y p) 

such that: 







"Xi" 






=U- 




+V 


-Z m _ 




_x p _ 





and 

Vk{1 pjyi^zimod N } 



[01 25] Since the from and to sets are the same and are finite, it is clearly necessary and sufficient to show that <pu 
is injective in order to show that q>u is bijective. 

[0126] Let Z denote the set of signed integers, then bijectivity of <p y is equivalent to the injectivity definition, i.e.: 
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V(q,,q J ...q p )e Z p 



V(xi, x 2 ... x p )e |Ni|x |N 2 x..,x N, 



v (yi, yi- - y P )e |Ni|x|n 2 |x..,x n. 





"xf 




"yi" 




q,-N,~ 




~xf 


"yr 






+v = u- 




+V + 




=> 






J 


_x p _ 




-yp- 










_y P _ 



*o [01 27] Then by a simple subtraction, by using linearity of U, and by making some variable change on x r yj it comes 
that bijectivrty of <pu is equivalent to 



15 



20 



V(q I) q2...q P )6Z' > 

V(xi, x 2 ... x p )€[-N,+l ,N,-l]x...x[-N p +l ,N P -1] 



"xi~ 




"q t Ni 


=> 






~0" 


_x D _ 








_Xp_ 




_0_ 



[0128] Where in the equation above all [a,b] denotes an interval of signed integers comprised between a and b 

inclusive. 

25 [01 29] Now let us assume that U is a sub-diagonal matrix, with a diagonal of 1 or -1 , that is to say: 



30 



v(ij)G{1 f ...,m} 2 ,j>i^u u = 0 

and 

IV i € {1 m}, u,,i = 1 oruii= -1 



(with of course all the elements of U integers). 

[0130] If such a condition is fullfilled then U is clearly inversible and the inverse U" 1 of U is also a sub-diagonal 
35 matrix with a diagonal of 1 or -1 and all elements still integers. This is trivially derived by the Gauss pivot algorithm 
applied to inverse matrix U. 
[0131] The equation 



40 



45 



u 







"qvNi" 


-Xp_ 




-q P Np. 



can be written: 



50 


"xf 




"qvNf 






= ir 1 . 






_Xp_ 




-q P Np_ 



55 



[01 32] Let us denote u'y the elements of U * 1 =[u' ,j] 

[0133] IT 1 is a subdiagonal matrix with all element integer and a diagonal of 1 or -1 . 
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• x 1 is clearly null, because x 1 = u , 1 M -q 1 • N 1 , q 1 and u' 1p1 are integer, and 0 is the only multiple N 1 in [-N^l , 
N r 1] 

5 • Now for all j with 1 sj<p, if vi 1^4 Xj = 0, then we have x R1 = u' j+1J+1 • q j+1 • N j+1 , and similarly we can deduce 
that Xj +1 is also null. 

[0135] Now let us assume that U*U (2) • U (1) , where U (1) is a subdiagonal matrix with all elements integers and a 
diagonal of 1 or -1 , and U (2) is a diagonal matrix whose for all i, i th diagonal element is primary with N|. 
10 [0136] Let us assume that for some (q 1( .... q p ) integers, and for some (x 1( ...Xp) in 
[-N t +1 ,N r 1]x...x[-N p +1.N p -1]wehave: 



15 



u 



"xr 




"qvNi" 


_x p _ 




_q p N p _ 



20 



[0137] Let us call 



25 







"X1 




= u (1) - 




-y P - 







30 



we have therefore: 



35 



U (2 > 







"qrNr 


_y p _ 







40 



[0138] That is to say: 



45 



r ( 2 ) i 
Ui.vyi 



qvNi 



50 



[0139] For all i from 1 to m, Nj divides uC 2 ) • yj, and is primary with u®} by definition of (J* 2 )' then, according to 
Gauss theorem, it also divides y s . In other word we can find some q \ , q^ integers such that: 



55 
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"yf 




"q;.N," 


-y p _ 




-q P Np. 



that is to say: 



10 



















»5 


-Xp. 







20 



[0140] Now is a sub-diagonal matrix of integers with diagonal elements all equal to 1 or to -1 . Then, as we have 
already shown it above, this involves that 



25 



Xi 



LXpJ 



= 0, 



and since fat q p ) were any integers, and (x 1( ...Xp) any elements of [-N ^ , N r 1]x...x[-N +1 , N -1] , this involves 

so that <pu is injective, and therefore bijective, thus concluding the proof. H 

ELEMENTARY FACTOR MERGE <p n 
[01 41 ] In order to get back to 



35 



s w =m 



40 



in the end, some x product of sets must also sometimes be merged. The elementary factor merge <p n function merges 
the two first factors of the first term of S( n_1 >. In fact the same principle could be used on any pair of factors of any term 
that have at least two factors, but this is not needed for the definition thanks to the term and factor permutation functions 
defined above. 
In that case: 



45 



50 



q n • q n . 1 , i.e. same number of terms 

Pi,n-1 * 2 and v , „ Pin=P 
vi€{2,... ( q n ; 



. ,,,, . i^r 1 - » one factor fewer in first term 
T, n =T, , i.e. terms other than the first one are unchanqed 



V J ^{3 p 1jH } N^- n^* ,; , i.e. the p 1>rK) -2 last factors of Tf-D are unaffected 

A and B defined by A^ N(";]) , B= Nft, 1 ) , and thus A • B= 
a bijection 



55 



h: 



0x[l 
(X1,X 2 ) h+ 



AB 



h((xi,x 2 ) ) 



<p n is built as: 
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<Pn: 



S <n-1) 



S (n) 



10 



15 



with 
x = (t,v) 
and 



v= ( x i;VV,) 



( 1 l( h(x 1 ,)C 2 ),X3._X pi/vi)) ift=1 

x otherwise 



[0142] The h bijections used by the elementary merge function can be defined as follows: 



20 



hjAl><fBl-> l a»B) main factor merge functions 



25 [0143] There are many possible functions that all are computable on any current processor. 
[0144] Here is one example: 
In the following for all x and y integer such that y>0: 



30 



35 



40 



x div y stands for the quotient of Euclidian division of x by y, 
x mod y stands for the remainder of Euclidian division of x by y. 



Note that vx vy > 0, (x mod y) e g], even for x ^ 0 
The DIV-MOD factor merge function: 



(X1,X 2 ) 



AB 



B*i+x 2 



[0145] This function is clearly bijective and its inverse function is the div-mod factor split g bijection exemplified in 
the definition of the elementary factor split <p n function: 

[01 46] Note: it is not needed to define a mod-mod factor merge h function that would be the inverse of the mod-mod 
45 factor split g function defined for the elementary factor split <p n> that is to say: 

g(x) = (x mod A, x mod B) . 

g" 1 ( (a,b)) = b + ( {(a-b) • B') mod A) • B 

50 

[0147] As a matter of fact the <p n that would result from this mod-mod factor merge h would be the composition of 
an affine-mod <p n of matrix defined by block : 



55 
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^ 



u (1) 



70 



15 



U= 



Id 



0 H 



O v 



B' 0 
0 1 



Id 



O h 



1 -1 
0 1 



20 

and of null vector V, and of the elementary factor merge function using the div-mod factor merge h function. 

[01 48] Above Id, Ov, and O h are respectively the identity (m-2)x(m-2) matrix (diagonal of 1 ), the null (m-2)x2 matrix 

and the null 2x(m-2) matrix. 

25 ELEMENTARY TE RM MERGE y D 

[0149] In order to get back to EO in the end, some * sum of sets must also sometimes be merged. The elementary 
term merge <p n function merges the two first factors respectively of the two first terms of S (n * 1 >, when the two first terms 
have the same subsequent factors, that can therefore be factorised. In fact the same principle could be used on any pair 
30 of factor of any pair of terms in provided that they can be factorised in the form 



(S*[b])x[a3x...x[S]. 



[0150] In that case: 

• q n -i^2 and q n = q n . n -1 . i.e. one fewer term 

vj e {3 p rv1 } t[" } = Tj n * , i.e. remaining terms unchanged 

40 • p 1M5 p. M sp 1ni i.e. the two first terms of S< n " 1) have the same number of factors, that is also that of the first 
term of 

• vj <= {2 p 1 ^ } N?j 1) = N£ n j 1) - that js t0 sa y tne to® first term can be factorised 

• A and B defined by' A 4 NO};}) , B 4 n^) , and thus A+B= N^ 

• a bijection k: 0 & [| IA+BI 
45 <p n is built as: 
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cp„: S (M > -+ S (n > 

s * h» fti.^t.x,)),^ x pw )if t=1 ort=2 

l(t-1,v) otherwise 

With 

w x = (t.V) 

and 

v= ( x v x 2- > Wt) 

[01 51 ] The term merge k bijections used by the elementary term merge function can be defined as follows: 
20 rt [A|^[S] -> [ A+Bj MAIN TERM MERGE FUNCTIONS 

[01 52] There are many possible lunctions that all are computable on any current processor. 
25 [01 53] As an example of k, the mod-threshold term merge function will be described hereafter : 

[01 54] In this function, a merge function is made that is the inverse of the mod-threshod term split function that was 
already explained. 

[01 55] Similarly we have two constants C and T such that: 
30 0< T < C 



35 



40 



(t,v) k> x defined by the algorithm below 



[01 56] Three constants x M , v M and t M are derived from A, B and C according to the following definitions and formu- 
las (the formulas are the same as for the mod-threshold term splitting function): 

x M is the greatest x while neither 0 nor [B] are empty. 
45 t M is the tag of the set (0 or [§)) not yet empty. 

v M is the greatest value of the set not yet empty while both sets are not yet empty. 
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A 

x M = max 



X € 



/7-T + c<A ifc<T 
A+B J and , where / - x div C. and c = x mod C 

V / (C-T) + c-T < B otherwise 



if((x M +l)mod C)>T 
otherwise 



10 



15 



f 



V 



max S 



max 



'/3x 



'/Bx^xm, 



6=xdivC 
and 

Vr =x mod C 

6=xdivC 
and 

\c =x mod C 



'c<T 
and 
= /-T + c 

'c>T 
and 

ly = /-(C-T) + c-T 



if t M = 1 



► if t M 38 2 
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[01 57] Then the algorithm that maps (t, v) onto x is the following: 

if v < v M or t*tM then 
begin 

remark Neither [a] nor [b] are yet empty, then discrimination is 

based on tag 
ift = 1 then 

x := (v div jyc + (v mod T); 
else 

x := T + (v div (C-T))C + (v mod (C-T)); 

end 
else 
begin 

remark [a| or \b\ is already empty, then no more discrimination, the 

remainder of the set not yet empty is mapped to |A+B 
x := v - v M + x M ; 

end 



55 



[0158] IMPLEMENTATION MATTERS 

[0159] In fact the method of the invention is not tied to any particular implementation scheme. The method only 
defines mathematically a permutation, with a definition allowing less tabulation at the expense of more processing. 



28 



09/10/2004, EAST Version: 1.4.1 



EP 1 039 646 A1 



[0160] As was already suggested above, the computation of I" 1 from its q>k«<Pk-i-.o<Pi definition might be done by a 
general purpose processor, typically during connection time, in order to prepare a look-up table to be used during inter- 
leaving. But as well, this computation might also be done by dedicated hardware circuitry, computing on the fly 
from its q>k°9k-i definition in order to interleave the (i+1 ) th symbol of data when this symbol is to be interleaved. Typ- 
5 ically a hardware circuitry could use a pipeline architecture where each step in the pipeline more or less correspond to 
one <p n , or one compound <p n+p o...o<p n . This kind of architecture is applicable because in order to compute on the fly r 1 (i), 
what matters is not the total computation time of I" 1 (i), but the rate at which they can be computed. This is because the 
T arguments of l" 1 (i) come in a predefined order. 

[0161] When a hardware architecture is used, the operations implemented are not necessarily using the mathemat- 
10 ical definition of <p n . In order to reduce the circuit size, there can be a mixture of arithmetic units, and of ROM tables. For 
instance a first x splitting step can be done arithmetically, then the intermediate steps can be done by look up table, and 
a final x merging step can be done arithmetically again. 

[0162] When a general purpose processor is used, neither are the operations implemented necessarily directly 
using the mathematical definition of q> n . There might be some optimisations : 
15 [0163] For instance when N is a power of 2, then splitting [FJ] into a product of sets whose respective element num- 
bers are also power of two does not need any division or modulo computation since it can be done by simple bit manip- 
ulation such as ORs, ANDs and bit shifts. For instance when mapping |256| to [8] x [32] by the div-mod factor split 
function, then the image of x is (y,z), where if X=b7b6b5b4b3b2b1b0 in binary notation, then y=b7b6b5 and 
Z=b4b3b2b1b0. To take it in more practical terms, if x=101 10001, then y=101 and z=10001. 
-20 [0164] Also still in the case of a number N of elements that is a power of two, and when there is only one term (or 
when only one term of size N in the & sum is considered), some factor permutations are implementable by classical 
processing units instructions such as bit rotation or quartet swap. Let us for instance consider for N=256 the 16x16 rec- 
tangular inter leaver. The 16x16 rectangular inter leaver comprises in its <p n compound definition a factor permutation 
that is swapping the two f}6] factors. In fact, the interleaver permutation can be shown to be not much more than a quar- 
ts tet swap which can be done in one instruction on some processing machines. 

NON-UNIQUENESS OF DEFINITION 

[0165] The inventive method provides means to define an interleaver using less tabulation at the expense of more 
30 processing, so that to reduce ROM size requirement. 

[01 66] The definition of the interleaver provided by the method might be not unique, since for instance: 

consecutive splitting steps can sometimes be done in several di ffere nt ways; for instance a first way to map |140| 
onto [7] x [4] x (5) is a first div-mod splitting function cp 1 mapping |140| onto GO x [^.followed by a second div-mod 
35 split 94°<P3«P2 mapping [20] onto [4]x[5] (q>2 and tp 4 are the factor permutation function respectively bringing [20] in 
first position, and then [7] in first position again); a second way is a first div-mod splitting function <p'-| mapping [go] 
onto ilx[5], followed by a second div-mod split <p' 2 mapping HH onto [3x0; both ways yield the same mapping 

(^•wi^Wili 

- the order of compound of two functions whose action is limited to different factor or terms can be inverted. 

40 - it is possible to build a function <p n and its inverse by the inventive method, then you can always insert a compound 
of several <p n into the compound defining l' 1 , provided that the compound of these several (p n sums up to the iden- 
tical function, and has no overall effect, for instance when r 1 =<p k o...«<p 2 °^i is an interleaver over El. and 
N=L • C , then if <p k+1 js defined as the div-mod factor split function that maps [R] to [U| x [c] and q> k+2 as the div- 
mod factor merge function that maps [Qx[c] to M, then I" 1 can also be defined as I" =<p k+2 °<Pk + i«<Pk-- ^2*^1 

45 because the compound q>k+2 o( Pk+i has no overall effect. 

the classification of the <p n might not be unique, some cp n can be found under two categories. For instance when all 
the N (f 1 ] are equal then a "factor permutation" <p n can also be defined as an "aff ine-mod" <p n with V=0 and U a per- 
mutation matrix. 

so EXAMPLE OF INVENTIVE METHOD APPLIED TO KNOWN INTERLEAVERS 

[0167] The inventive method allows to build interleavers by mixing tabulation (or processing not defined in the 
invention) and processing (defined in the invention). Here is shown that some well known interleavers can be defined 
with the inventive method to be implemented in an interleaver device according to the invention. 

55 

1 - THE RECTANGULAR INTERLEAVER 

[01 68] According to the invention, the classical rectangular interleaver I of depth N=L • C with L lines and C column 
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block is written as: 

,1 =9 3 >2>1 

5 Where: 

• <pi is the div-mod elementary factor split function splitting El to [Qx[c] 

• <p 2 is the factor permutation function for o^(2,1), that is to say cp 2 is mapping on [C]x[L] 

9 3 is the div-mod elementary factor merge function merging \C\ x [Q to El The rectangular interleaver is illustrated 
10 on figure 6. 

2 -THE DIAGONAL INTERLEAVER 

[01 69] The usual diagonal interleaver I depth N=L • C , with L lines and C column block, where lines are written first 
15 and then diagonal are read, when j th diagonal begins at the first element of j* 1 line is written as: 

r 1 = ( p 3 o (p 2 o (p 1 

Where: 

20 

<Pi is the div-mod elementary factor split function splitting [HI to [C]x[Cl 
9 2 is * ne affine-mod function with V=0 and 



30 • <p 3 is the div-mod elementary factor merge function merging [Qx[c] to El 
3 - THE TRUNCATED RECTANGULAR INTERLEAVER 

[0170] Truncated rectangular interleavers are used in MIL interleavers which are defined in ETSI SMG2 UMTS-L1 
35 EG Tdoc 98/273 "Multi-stage interleaving (MIL) method for Turbo codes" by "NTT DoCoMo". In a truncated rectangular 

interleaver the N input data is written along lines in a block of L' lines and C* columns, with C* which is not a divider of 

N (and L' great enough for L' • C £ N). Then the data is read out along columns. 

[0171] The truncated rectangular interleaver is illustrated on figure 7 with N = 20 and C = 6. 

[0172] "Rie truncated rectangular interleaver is defined by the method of the invention as the sum of two rectan- 
40 gular interleavers. 

I " 1 =<P 8 °(<P 7> 6>5)°( ( l ) 4>3 0( P 2)°<P 1 

Where : 

45 

9-I is the mod-threshold elementary term split function, where C = C* and T = N mod C . 
91 El to[N]* [N^J with 

N 1 = (N div C+1) • (N mod C), and 

50 

N 2 = (N div CV(C*-(N mod C 1 )) 
(<P4°<P3>2) j s the first rectangular interleaver that is applied to [ffj] with parameters: 
55 L = (NdivC , +1),and 

C = (N mod C) 
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(<P7°<P6 0( P5) » s the second rectangular interleaver that is applied to [FJ^] with parameters: 

L=(N div C*), and 

5 C = C - (N mod C) 

(p 8 is the term merge function with C = N .,+N g and T = N A . 
4 -THE MIL INTERLEAVER 

10 

[0173] The MIL interleaves defined in ETSI SMG2 UMTS-L1 EG T doc 98/273 "Multi-stage interleaving (MIL) 
method for Turbo codes" by "NTT DoCoMo" can also be constructed by the method of the invention. The MIL interleaver 
is defined recursively, and the method of the invention is by essence recursive. 
[0174] A MIL interleaver of depth N = L * C is defined as the succession of the following steps: 

15 

1) First write the input data along the lines of the rectangular block of L lines and C columns; 

2) Then interleave each line by either a rectangular interleaver, or a trunctated rectangular interlaver, or a MIL inter- 
leaver; and 

3) Finally read the block by columns. 

20 

[0175] The first step consists of (9 2 °<Pi). where <p 1 is the div-mod factor split function applied to El and mapping it 
to [E]x[c], and 9 2 »s the distribution function that maps [EJxED to 



25 




L times 



30 where each term consists of a line of the block. 

[0176] The second steps consists in applying a rectangular, a truncated rectangular or a MIL interleaver to each 
term. As already seen, rectangular and truncated rectangular interleaver can be generated by the method of the inven- 
tion, then this second step is also within the invention. 

[0177] The third step consisting of reading the block by column can be simply written as (<Pk°<Pk-i°SPk-2). where <p k _ 2 
35 is the factorisation function mapping 




L times 

40 

to \Qx\c\, (p k _-, a factor permutation mapping to [C]x[L], and <p k the div-mod factor merge. 

[0178] In conclusion MIL interleavers are a subset of the interleavers that can be generated by the inventive 

method. 

45 

5 - OTHER INTERLEAVERS 

[01 79] The inventive method is open to embed other interleavers into the inventive interleavers by means of the ele- 
mentary interleaving function. After decomposing 

50 

S<°>=[n] 

into a sum of products of sets of the 0 form, an interleaver not defined by this document can be applied to only one of 
55 this factor set : this is what is called elementary interleaving in this document. 

[0180] In the definition of many existing interleaver a similar method is also used, for instance in ETSI SMG2 
UMTS-L1 EG Tdoc 98/337 "A General Turbo Interleaver Design Technique with Near Optimal Performance " by " 
Hughes Network Systems", we can read that such a process is used : 
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• writing the symbol in a rectangular block along the rows: this corresponds to a div-mod factor split of El to [L]x[C] 
then operating on each row separately: this corresponds to a distribution of [Qx[C] product in a sum of L terms 
* ... \C\, each term being a row, 

shuffling each row by some rule derived from Galois field arithmetic: this corresponds to an elementary interleaving 
5 applied to the term of which the row is consisting. 

permutating the rows : this can be expressed either as a term permutation, or as factorising to [Qx[c] again, and 
making an elementary interleaving bearing on the [Q dimension. 

reading by column: this is also interpreted in the inventive method, we first swap the two factors in [L]x[C] by a factor 
permutation, then make a div-mod factor merge. 

10 

[0181] The invention brings some methodology allowing to define an interleaver in a unified formal language. This 
would ease specification of interleaves by defining the compound or sequence of the <p n . Such a specification would 
be without ambiguity, and even possibly in a language that is directly machine processable. Such a sequence can be 
coded in a dedicated binary format, for instance a TLV format, where T (tag) gives the type of (p n (factor split, term split, 

15 affine-mod, etc..) and the LV (length+value) gives the characteristics of <p n , for instance if <p n is affine-mod, the matrix 
U and the vector V. Such a coded format can be used to defined a great variety of interleaves in a device at a very low 
ROM memory cost. Also the same format can be used on an Interface to negotiate which interleaver is to be used. 
[0182] The term split is also of a great interest. Many arithmetic inter! eavers (e.g. GF interleaver as defined in ETSI 
SMG2 UMTS-L1 EG Tdoc 98/337 "A General Turbo Interleaver Design Technique with Near Optimal Performance " by" 

20 Hughes Network Systems", for N=2 m ) are primarily defined for some N that has good properties. When an interleaver 
with an N that has not such good properties needs to be built, then an arithmetic interleaver with N*>N and N' having 
the good properties is built, and then punctured down to N. Puncturing an interleaver is not done at a very high process- 
ing power when it consists in modifying a RAM table. However these forces an implementation where first the inter- 
leaver is prepared in a RAM table, and then this RAM table is used. Implementation where the interleaver function is 

25 computed real time by an ASICs cannot be realistically considered with puncturing. Thanks to the term splitting function 
instead of searching for N' with good arithmetic properties and N'>N, one can find N' with good arithmetic properties 
such that N'<N, and then apply the arithmetic interleaver to a term of size N', thus avoiding any puncturing. 
[0183] Puncturing of an interleaver is shown on figure 8. An interleaver l 8 of depth 8 is punctured down to a depth 
of 5 to form an interleaver l 5 . In the upper part of the figure a table where the values of l 8 " 1 are shown. The table entry 

30 number (above the table) is the value of i. and the entry value (inside the table) is the value of l 8 " 1 (i). In other words we 
have l 8 * 1 (0)=3, l8 _1 (1)=0, etc.... Then all the entries such that l 8 " 1 (i)^5 are to be punctured. These entries are marked 
on the figure by a triangle. Puncturing consists in moving some part of the table, for instance the entry in position 4 is 
moved to position 2 so that l 5 (2) = l 8 (4) = 1 . 

[0184] After puncturing, the 5 first entries of the table contain only numbers ranging from 0 to 4. 

35 

Claims 

1 . Interleaver device for interleaving a data set comprising: 

40 - processing unit (10) including a data processor (16) for running an interleaver (I" 1 ) under the control of driving 
means (18); 

- input means (12) for inputting the data set to be interleaved and output means (14) for outputting the inter- 
leaved data set; 
characterised in that: 
45 said driving means (1 8) include: 

mapping processing means (20) for performing a set of bijective elementary functions (<p n ) and supplying 
a mapping of the interleaver to the data processor (1 6) for interleaving the data set according to this map- 
ping; and 

so j n terleaver definition means (24) for supplying said mapping processing means (20), with a definition of 

said interleaver (I" 1 ), expressed as a compound function (9^... o^-j) of elementary functions (<p n ), each 
elementary functions coming from said set of bijective elementary functions (<p n ), for said mapping 
processing means (20) to perform each of said functions, compounded according to the interleaver defini- 
tion (I' 1 ), and so providing to the data processor (16) said mapping according to which the data processor 

55 (16) interleaves the data set and supply the interleaved data set to the output means (14). 

2. Interleaver device according to claim 1 , characterised in that the driving means (18) includes a read only memory 
(22) from which the set of the bijective elementary functions (<p n ) can be derived. 
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Interleaver device according to claim 1 or 2, characterised in that the interleaver definition means (24) comprise 
means for receiving the interleaver definition (I' 1 ). 



Interleaver device according to any one of the preceding claims, characterised in that said set of elementary func- 
tions (<p n ) includes at least the following functions whose origin and destination sets can both be expressed as a * 
sum of x products of factor sets that all are sets of successive integers: 

an elementary term split function, which is splitting, into a * sum of two sets, a factor of a term of a set; 
an elementary factor split function, which is splitting, into a x product of two sets, a factor of a term of a set; 
an elementary permuting function, which is permuting the elements of a factor of a term of a set; 
an elementary aff ine-mod function, which is acting on a term of a set as if it was a subset of a vector space 
whose dimension is the number of factors, by an aff ine mapping, followed by taking the modulus of each co- 
ordinate by the size of the respective factor set; 

an elementary factor merge function, which is merging a x product of sets, into a set by merging two factors of 
a term; 

an elementary term merge function, which is merging a ^ sum of sets into a set by merging two factors 
respectively of two terms of a set, when the two terms have the same other factors. 

Interleaver device according to claim 4, characterised in that said set of elementary functions further includes 
the following functions whose origin and destination sets can both be expressed as a * sum of x product of factor 
sets that all are sets of successive integers: 



a elementary term permutation function ; and 
a elementary factor permutation function; 
and in that: 



- the elementary term split function is splitting, into a ^ sum of two sets, the first factor of the first term of a 
set; 

the elementary factor split function is splitting, into a x product of two sets, the first factor of the first term 
of a set; 

the elementary permuting function is permuting the elements of the first factor of the first term of a set; 
the elementary factor merge function is merging a x product of sets into a set, by merging the two first fac- 
tors of the first term of a set; 

the elementary term merge function, which is merging a ^ sum of sets a sets into a set, by merging the 
two first factors respectively of the two first terms of a set, when the two first terms have the same subse- 
quent factors. 

Interleaver device according to claims 4 or 5, characterised in that the elementary factor split function is splitting 
one factor |A* B| of a term in the origin set by the following mapping: 



x ^ (xdiv B, x mod B) 



Interleaver devic e according to claims 4 or 5, characterised in that the elementary factor split function is splitting 
one factor |A* B| of a term in the origin set by the following mapping: 

|Xb] -+ [a]xU 

x (x mod A f x mod B) 



Interleaver device according to anyone of claims 4 to 7, characterised in that the elementary factor merge function 
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is merging two factors 0 and H of a term in the origin set by the following mapping: 



0*1 
(X1,X 2 ) 



A-B 



Bxi+x 2 



10 9. Interleaver device according to any one of claims 4 to 8, characterised in that the elementary affine-mod function 
is defined by the affine function of matrix U which is the product U ■ U (2) • U (1) of two matrices U (1) =[u ■ ^] , and 
U (2) =[u^] such that 



15 



20 



25 



30 



U (1) is a sub-diagonal matrix with a diagonal of 1 or -1 , that is to say: 



V(i,j)e {1 m} 2 j>i=*u ( ,y=0 

and 



Vk{1 m}, u ( J' = 1 or uYj' = -1 



.w - 



is a diagonal matrix with for all j, j* 1 diagonal element primary with N <A> that is to say: 



V (r,c)€ {1,...,m} 2 r*c => u ( r 2 c= 0 
and 

V je {1 f ... t m>, u'jj and U\f greatest common divider is 1 



1 0. Interleaver device according to any one of claims 4 to 9, characterised in that the elementary term split function is 
35 splitting one factor IA+BI into |A| * H by selecting the tag of the image based on a com pariso n to a constant thresh- 
old of the modulus of its IA+BI input by a some fixed constant, at least for some of its IA+BI input values. 

1 1 . Interleaver device according to any one of claims 4 to 1 0, characterised in that the elementary term merge function 
is merging into |A+B| two factors 0 and @ taken in terms with same respective other factors, such that its inverse 

40 bijection selects the tag of the image in 0 * |g based on a comparison to a constant threshold of the modulus of 
its IA+BI input by a some fixed constant, at least for some of its IA+BI input values. 

12. Method for interleaving a data set according to an interleaver (I 1 ), characterised in that the interleaver (l~ 1 ) is 
defined as a compound function (q> k o...o<p 1 ) of bijective elementary functions (<Pn), each elementary function (tpj 

45 coming from a set of bijective elementary functions (<p n ), and in that it includes the steps in which the elementary 
functions, compounded according to the interleaver definition (I" 1 ), are subsequently applied to the data set for pro- 
viding the interleaved set of data. 
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